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Abstract. - We revisit the problem of compaction of a column of granular matter exposed to 
discrete taps. We accurately control the vertical motion of the column, which allows us to vary 
the duration T and the amplitude A of single-cycle sinusoidal taps independently. We find that 
the density of the material at the reversible branch depends both on A and T. By comparing 
the densities on the reversible branches obtained for a range of values of T, we find that we can 
collapse all data when plotted as function of A/T, which scales similar to both the liftoff velocity 
and the time of flight of the packing. We further show that switching between states obtained for 
different A and T, but chosen such that their densities on the reversible branches match, does not 
lead to appreciable hysteresis. We conclude that the appropriate control parameter for sinusoidal 
tapping is not the peak acceleration F A/T^, as is usually assumed, but rather FT ^ A/T. 
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Everyday experience tells us that granular materials 
ompact when vibrated — think of tapping a tin of coffee 
powder in a can. Early experiments in the 1950s indicated 
that tapping can induce both dilation and compaction of 
a packing of granular materials, depending on the details 
of the tap [1]. More than a decade ago, a series of exper- 
iments probed the compaction of dry granular materials 
(glass beads) in a narrow tall tube which was tapped verti- 
cally [2-4] . Starting from a loose packing, the density was 
observed to exhibit slow, logarithmic growth as a function 
of the number of applied taps. Memory effects, in which 
the density evolution depended on the tapping history, 
were also found [5]. 

Eventually however, in all experiments a state was 
reached where the packing density depends only on the 
tapping strength and not on the history [6]. From such 
a state the so-called reversible branch could then be ob- 
tained repeatedly and reversibly by slowly ramping up and 
down the tapping intensity. The transient phenomena and 
memory effects all occur along the so called irreversible 
branch, along which the system would evolve before reach- 
ing the reversible branch. 

In later studies of compaction in much wider contain- 
ers, convection was found to be important [7-9]. In these 
experiments the temporal evolution of the density on the 
irreversible branch was found to be different [6-9], but 



again the same reversible branch was found. 

Here we address the following question: What is the 
appropriate control parameter that characterizes the taps? 
A widely used characterization of taps is the ratio of their 
peak acceleration and the gravitational acceleration, F. 
Certainly the peak acceleration is important in that it 
allows to distinguish between tap strengths where no liftoff 
of the packing occurs, for F < F* ^ 1, and taps where that 
does happen [7] . Recent numerical work on the irreversible 
branch dynamics of compaction has suggested that the 
dimensionless acceleration parameter is not appropriate 
for rescaling the data [10]. Moreover, recent numerics [11] 
of compaction under sinusoidal driving vibrations indicate 
that the vibration frequency also influences the reversible 
branch, and similarly, supporting evidence for the role of 
tap duration can already be found in the observations of 
Macrae et al. [1]. 

Here we address the question of the appropriate control 
parameter experimentally, by studying the packing den- 
sity on the reversible branch in experiments in which we 
expose granular packings to single cycle sinusoidal discrete 
taps with r > 1, where we control and vary both the tap 
amplitude A and its duration T (Fig. [T]). For sinusoidal 
taps as used here, F can be given in terms of A and T: 
F = Auj'^ /g^ with uo = 27r/T the radial frequency, and g 
the gravitational acceleration. Since we have precise con- 
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Fig. 1: (a) Sketch of the experimental setup, (b) Sketch of 
the waveforms of z{t), z{t) and z(t), iUustrating the offset and 
hnear slopes present in z{t), stemming from continuity require- 
ments, (c) Comparison between desired waveform (grey) and 
actual waveform (black) of the vertical acceleration z{t). (d) 
A typical measured position (black) and acceleration (grey) 
signal. The arrow indicates the phase at which the packing 
fraction is measured. 



granular media subjected to vibrations [12]. An 0.5 me- 
ter long and flexible aluminum rod couples the tube and 
the shaker (see Fig. [l^). This rod eliminates the neces- 
sity of excessively accurate alignment of the axes of the 
air bearing and the shaker. 

The vertical acceleration of the tube is measured with 
an accuracy of 10~^g, by a combination of two accelerome- 
ters (Dytran E3120AK and a modified ADXL320EB). For 
consistency checks, we occasionally monitor the position 
z{t) of the glass tube with an inductive position sensor 
(Messotron WLH50, 10~^m resolution) — the accelera- 
tion measurements are far more sensitive and these will 
be used in the feedback scheme described below. 

Tap — In order to create a tap where both the du- 
ration and strength can be controlled independently, we 
need to determine the waveform z{t). We denote the du- 
ration of the tap by T = 27t/uj^ the period of the tap by 
P (we will fix it later at 1 second), and have the tap start 
at t = 0. We demand than that z(t), z{t) and z{t) are 
continuous, that J dtz{t) = and J dtz{t) = 0, and that 
z{0) = 0. These requirements severely constrain the tap 
— for example, taking for z{t) a single sine cycle followed 
by a period where z = makes z{t) discontinuous. Start- 
ing from a single sine cycle for z{t)^ obtaining z{t) and 
z{t) by integration, and fixing the integration constants 
so that the continuity conditions are full filled, we find the 
waveform summarized in Table [l] (see Fig. ^p) : 



trol over the acceleration signal, we can vary both F and 
uj independently (see Fig. [T]3). 

For given T, we obtain similar reversible branches 
(/)rev(r,^) as were observed before for fixed T, both for 
a bidisperse glass bead mixture and a bronze powder, and 
the reversible densities 0rev(r,T) depend both on T and 
F. We find that all data can be collapsed in good approxi- 
mation by plotting the packing fraction as function of FT, 
which is close to the liftoff velocity, or similarly, the time 
of flight of the granular packing. In addition we probe 
what happens when we alternate the taps between pairs 
of different F and T for which ^rev(r,T) are equal, and 
we find no appreciable hysteresis. 

The evidence that we will detail below will show that F 
is not the appropriate control parameter for compaction. 
Our results suggest that the density of the reversible 
branch is controlled by the product of F and large T. 
Hence, by maximizing F T instead of F, a wider range 
of agitation strengths can be reached. 

Experimental Methods. — Setup. We study granu- 
lar compaction in a glass tube (diameter 26 mm, height 20 
cm) filled with grains (typical filling height 10 cm). The 
glass tube is shaken vertically with a shaker (VGIOO, VTS 
systems), driven by a commercial audio amplifier (Crown 
CEIOOO). The vertical motion of the tube is guided by an 
air bearing (0 1", New Way) which is levelled. Levelling 
eliminates heaping, the unwanted tilting of the surface of 
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Table 1: Waveform of the tap. Here the integration constants 
equal bi = A(jo{P - T) / P and 62 = auo T/P. 



In order to create such a tap, a feedback algorithm is 
used to adapt the Fourier components of the signal fed into 
the audio amplifier such that the measured acceleration 
signal converges to the desired waveform. This procedure 
is carried out once, before the start of the experiment, and 
the shape of the output wave is not changed during the 
actual experiment; it is only multiplied by a scale factor 
to set its overall amplitude. We have checked that the 
system is sufficiently linear so that re-calibration of the 
waveform is not necessary when F is varied. 

Taps are applied once every second, so P = 1; this 
allows the packing to come completely to rest before the 
start of a new tap. In Fig.jl]^, the desired waveform for z{t) 
is compared to a typical measured acceleration signal — 
the waveform produced by the feedback scheme is in good 
agreement with the desired tap. In Fig. [l]i, a measured 
waveform for z{t) and accompanying z{t) are shown. The 
slow downward motion of the tube anticipated in Table [l] 
is clearly visible. This downward motion is not at constant 
speed due to low-frequency limitations in the electronics 
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Fig. 2: (Color online) The granular material used in the exper- 
iment (scale = 1mm): (a) glass beads, (b) bronze powder. Illu- 
mination direction: (c) for glass, from the back, (d) for bronze, 
from the top. Typical camera image (e) glass, (f) bronze. For 
bronze powder two cameras are used, image is typical for both. 

that drive the shaker. The accelerations associated with 
this non-constant speed are however always below so 
they do not influence the packing fraction measurably. 

Parameter Range and Grain Dynamics. As was estab- 
lished in [7], the dynamics of compaction dramatically 
slows down when the tapping peak acceleration F falls 
below a critical value F* which is typically slightly above 
1. This transition is accompanied by a change from a 
regime where the grains experience lift off for F > F* to a 
regime where they do not. In our experiments, the values 
of F range from 1 to 15, and by high speed imaging we 
have established a clear lift off and expansion of the whole 
column of grains for F > 2 — the precise value of F* in 
our experiments is likely somewhat smaller. The vast ma- 
jority of our data is therefore taken in the regime where 
the grains loose contact with the bottom plate and their 
neighbors. 

Packing Density and Material Used. The packing 
density is determined from the height of the granular 
column, and this height is measured with a camera that 
observes the packing from the side, as in Refs. [13, 14]. 
The camera is triggered for strobed image acquisition, 
such that the camera takes a picture just before the start 
of the tap: the trigger moment is indicated by the arrow 
in Fig. [l]i. A ring-shaped dark marker, whose height 
from the bottom of the tube is known (typically 13 cm), 
is attached to the tube, slightly above the maximum 
height of the bead pack. Depending on the grains used 
we use lighting from behind (for glass beads), or from the 
top (for bronze powder) — see Fig. [2] Determining the 
height of the packing in both cases amounts to counting 
pixels with intensity above a certain threshold value. 
The threshold value is determined from the histogram of 
a typical image. We verified that the imaging method 
applied gives a linear relationship between the amount of 



grains in the tube and the number of pixels in the gap. 

We employ two types of granular matter: a bidisperse 
mixture of glass beads (Pneumix, 1.6 and 2.3 mm in 1:1 
volume mixture, Fig.|2^), or monodisperse bronze powder 
(Acupowder, grade '12HP', ~ 1mm diameter. Fig. ^p)- 
Due to their different opacity, optimal lighting is differ- 
ent for these materials (Fig. |2]3-d). One more difference 
is that glass bead packings generally have a flat surface, 
while the bronze beads grow a small radially symmetric 
heap at their surface (this heap is likely due to convection) . 
One important difference is that the glass beads tend to 
be more sensitive to triboelectric effects, which result in 
overall, run-to-run variations of the absolute density (the 
trends in the variation of ^rev with F are not affected by 
such variations). We observe that such variations are ab- 
sent for the bronze powder. Importantly, our essential 
findings (control parameter ~ FT and no hysteresis) are 
seen with both materials. 

Our measurement method yields a high resolution; the 
noise in the determination of the height in the packing due 
to camera pixel noise alone is 30/im, this translates into a 
variation of (/) of ~ 0.03%. However, this method gives a 
relatively poor accuracy in the absolute value of (j). This is 
due to the cumulative effect of the errors in the determina- 
tion of the height of the marker, the inner diameter of the 
tube, the density of the beads, image calibration and the 
thresholding. The packing fractions stated in this paper 
therefore have an estimated systematic error of 1.2%. 

Experimental Protocol. The bottom of the tube is 
perforated and connected to a flow box so dry air can 
be pumped through the bead pack. We always start an 
experiment by fluidizing the packing with several pulses 
of dried compressed air. This procedure creates an initial 
packing density of order 0.60 ± 0.01. The number of flow 
pulses is generally not the same for each experiment; 
several pulses are applied initially to ensure a proper 
erasure of memory effects [2,14], after which pulses are 
applied until a packing with a flat surface is obtained. 
Airflow is turned off during all compaction experiments, 
which simply consist of observing how (p evolves while 
taps are applied. 

Transients & Steady State. — Fig. [3^ shows how 
the packing typically changes with the number of taps 
applied to the system, in an experiment where a packing 
of glass beads is subjected to 250000 taps with T = 33 
ms, and where the amplitude is altered between F ~ 2 
and F 6 after each 25000 taps. The packing fraction is 
only measured every 10 taps. There are several distinct 
features in this graph: (i) Initial transient: Initially the 
packing fraction evolves slowly towards a steady state 
value, nonmonotonically in the case depicted in Fig. [sJd. 
(a) Steady state: After the transient the system is able 
to reach a steady state packing fraction. The correspon- 
dence between (j) and F is also reproducible: changing tap 
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Fig. 3: (a) Evolution of the packing fraction as function of the 
number of taps. Here, we fix T = 33 ms, and F is alternated be- 
tween F = 2 and F = 6 every 25000 taps. Our data shows that 
after a short transient, the density can be reversibly changed 
from one value to another, (b) The initial transient compaction 
behavior on a logarithmic scale, for the range indicated by the 
dark grey bar in panel (a), (c) A typical compaction process 
between two steady states as occurs at t = 10^ as indicated by 
the light grey bar in panel (a). 
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Fig. 4: Reversible and irreversible branch for the glass bead 
mixture (See text). The main panel shows the reversible branch 
(/)rev(F,T) for different tap lengths T (+ 50ms, o 33ms, A 
16ms, * 14ms, □ 13ms, x 10ms). Clearly F alone is not suffi- 
cient to characterize the tapping. The inset shows (/)rev(F,T) 
for T= 14.3ms, and four sweeps in F (see text). Greyscale of 
symbols indicates measurement time; lighter points are later 
measurements. After an initial transient (black crosses), the 
densities become history independent. 



amplitudes repeatedly back and forth shows that for a 
particular tap amplitude the same steady state packing 
fraction is always obtained (Fig. |3^). (in) Amplitude step 
transient: When the amplitude of the tap changes, the 
packing fraction changes too, but not instantaneously. 
This time dependent process is usually not symmetric 
in the step direction: an increase of the tap amplitude 
usually almost immediately dilates the packing to the 
packing fraction appropriate for the new amplitude. 
A step down in tap amplitude requires compaction of 
the packing to reach the steady state packing fraction 
belonging to the lower tap amplitude. This process is 
usually slower. As can be seen in Fig. [Sj^, the increase 
in density from roughly 0.61 to the plateau value ~ 0.63 
takes of the order of 1000 taps, comparable to the initial 
transient duration. 

The initial transient is somewhat different in each ex- 
periment, and run-to-run fluctuations of details of ^(t) are 
considerable. For example, non-monotonicity seen in 0(t) 
in Fig. |3] is not always observed. This may either be due 
to the fact that the preparation of the packing in each ex- 
periment is slightly different, or to inherent fluctuations of 
the density evolution. We identify the transient behavior 
with the irreversible branch, and the steady state with the 
reversible branch. 

The strength of the fluctuations in the early evolution 
in (j){t) hinder a precise comparison of our results to the 
results from the Chicago and Rennes groups [2,3,9]. How- 
ever, our values for the densities on the reversible branch, 
^rev(r), are far more robust, and in the rest of the paper 
we focus on 0rev(r,T) in the steady state. 



Steady-State Density as a Function of F and 

T. — As we will show, our experimental protocol allows 
the determination of 0rev(r,T). At fixed T, this enables 
use to reproduce the reversible and irreversible branch 
in a 0rev(r, T)-plot similar to those found before in the 
Chicago and Rennes experiments by sweeping F up and 
down several times. 

In the inset of Fig. |4] we show that the packing density 
becomes well-defined on the reversible branch. In this 
experiment, starting from a low density packing, the tap 
amplitude is sweeped up and down four times, and the 
number of taps at each different tap amplitude is 4000. 
Every 20 taps the packing fraction is measured and each 
data point corresponds to an average over the resulting 
200 measurements. The irreversible branch is visible as 
the initial increase of (j) with F. After about 12,000 taps 
the reversible branch is reached: this branch is clearly 
shown in the four amplitude sweeps that all follow the 
same ^rev(r,T) relation, with 0rev decreasing with F. 

Effect of the tap duration. We will now explore the re- 
versible branch for a range of values of T (Fig. [4|. The 
reversible branches obtained in a series of experiments in 
which F is sweeped for a range of values of T are shown 
in the main panel of Fig. [4j The reversible branches for 
different tap lengths T clearly have the same overall form, 
but do not coincide — F alone is not the parameter that 
governs compaction. The spacing of the data and the fact 
that the functional form of ^rev(r,T) is similar for dif- 
ferent T, strongly suggests the data sets can be collapsed 
onto a master curve. Fig. Is] shows that ^ref(r,T) can be 
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Fig. 5: Data for the reversible branches for the glass bead Fig. 6: (/)rev(rT) for the rough bronze powder. Symbols: + 
mixture as shown in Fig.[4]can be collapsed well when plotted ^O^s, * 33ms, o 22ms, o 17ms, A 13ms , □ 10ms. Inset: Same 



as function of FT. 



data, 



,(F,T) for the different T. 



collapsed reasonably well by plotting ^rev as a function of 
the product FT. 

Bronze powder. For the bronze powder, we measured 
the reversible branch 0rev(r,T) for different values of the 
tap length T; the results are shown in the inset of Fig. |6] 
Note that the range of packing fractions obtained is far 
larger for the bronze powder, which may be attributed to 
the roughness of the particles [15]. Similarly to the glass 
beads, (p^ev collapses onto a master curve when plotted as 
function of FT — see Fig. [6j This shows that the details 
of the granular material used are entirely insignificant. 

Absence of hysteresis. In Ref. [5] memory effects were 
observed in the evolution of the packing fraction on the 
irreversible branch: two different initial conditions were 
prepared at a certain fixed (/), by compacting low density 
systems at different tapping strength. The time evolution 
of subsequent compaction differed, despite the fact that 
the initial packing fraction and tap strength were equal — 
hence not only the packing fraction, but also the history 
is important for the evolution at the irreversible branch. 

While such memory effects have not been seen, per def- 
inition, on the reversible branch, one may wonder if some 
more subtle hysteretic effects could arise there. In partic- 
ular, if the state of the system on the reversible branch is 
not fully specified by density, it might be that using two 
control parameters, subtle hysteretic effects not seen when 
only sweeping the single parameter F become apparent. 

Since data for ^rev as a function of F and T can be col- 
lapsed on a single master curve, it follows that F and T can 
be varied at the same time in such a way that ^rev stays 
constant. As we will show, we do not observe any appre- 
ciable hysteretic effects when switching between different 
pairs of (F,T) adjusted so that (/)rev(r,T) is constant. 

We measure the packing fraction of a bronze powder 
packing while we expose the granular packing to taps with 
a sequence of different T: 17,33,17,11 ms (Fig. [T^) and 
different F. This sequence is repeated three times; each 



F, T-pair was used for 10000 taps, the total number of 
taps applied in the experiment is 130000. The amplitude 
F (4.84, 2.88, 4.84, 7.32) for each different T, shown in 
Fig. [7]3, was fine tuned such that the resulting packing 
fraction stayed the same during the whole experiment; the 
product FT is constant to within 10%. This 10% varia- 
tion suggests that FT is likely an effective approximation 
of the ultimate order parameter that sets the density on 
the reversible branch — see below. The packing fraction 
evolution during the experiment is shown in Fig. [T]^. No 
appreciable transients are observed when switching F and 
T: there is no evidence for any appreciable hysteresis. Our 
data is further evidence that memory effects, which are a 
characteristic of the irreversible branch, do not play a role 
on the reversible branch. 

Interpretation. — How can we understand the rel- 
evance of FT in determining (j)rev^ We note that the 
liftoff velocity zi , the time-of- flight r/ and impact velocity 
Zi are all approximately, although not precisely, propor- 
tional to FT. Below we briefly discuss these parameters, 
but note already here that the experimental scatter pre- 
vents us from unambiguously determining whether these 
provide better scaling collapse than the simple FT scaling 
employed above. 

( i ) Liftoff velocity — The liftoff velocity zi is determined 
by calculating the velocity at the time where the accelera- 
tion falls below a threshold — ^F*. In case the packing ex- 
periences no friction with the container walls, this thresh- 
old equals the gravitational acceleration (F* = 1), but in 
general F* is somewhat larger, typically around 1.2. For 
given F* and waveforms as described above, one finds that 
in the parameter regime where liftoff occurs (F > F*), the 
liftoff velocity is given by 
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For most parameter values, the term within square brack- 
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when switching between different driving parameters that 
correspond to the same reversible-branch-density. 

The time of flight, the liftoff velocity and the impact 
velocity all scale similar to FT, which makes it impossi- 
ble to experimentally determine which parameter leads to 
the best data collapse. In experiments focussing on the 
transients in granular compaction [10], or on the hydrody- 
namic phases in vibrofluidize granular materials [17], the 
order parameter found was always proportional to the in- 
jected energy per vibration cycle. We suggest that our 
scaling supports the view that the energy injected is also 
the driving mechanism in compaction experiments. 



Fig. 7: T (a) and F (b) used in the constant packing frac- 
tion sweep. The packing fraction evolution in time is shown 
in (c); the fluctuations visible in (c) are of the order of the 
experimental noise. For details see text. 



ets is close to two — typically, T/P is of order of a few 
percent, and for F > 5, the square root term is within a 
few percent equal to one. The main deviations between 
the scaling of FT and zi occur for small F. Our scatter 
is relatively large in this regime, and in a plot of (/)rev as 
function of i^, the quality of the collapse is very similar to 
when FT is used. 

(ii) Time of flight — In rough approximation, the time 
of flight Tf is simply proportional to the liftoff velocity: 
Tf = 2zi/g. We have also calculated the time of flight 
numerically, assuming either the analytic expressions for 
z{t) and z{t), or using the experimental data available for 
z{t) and z{t), and taking drag forces into account (these 
affect the total time of flight of the packing [16]). Despite 
the presence of free fit parameters (take-off acceleration, 
drag coefficient), the resulting data collapse when Tf is 
used as scaling parameter is only marginal better than for 
FT or zi. 

(Hi) Impact velocity — The amount of energy dissipated 
once the packing comes to rest in the container is set by 
the impact velocity, and this could also be an appropriate 
choice of rescaling parameter. However, in the range of 
parameters explored, the impact velocity scales very sim- 
ilar to the time of flight, and the quality of the collapse 
does not improve. 

We hence conclude that while impact velocity or time 
of flight might give a marginal improvement in the data 
collapse, the much simpler rescaling parameter FT works 
essentially equally well. 

Conclusions. — By precisely controlling the shape of 
the taps, we find that the compaction of vibrated granular 
media is not only governed by F, the dimensionless peak 
acceleration, but that the tap duration plays an equally 
important role. We observe a collapse of the reversible 
branch with FT, for two different types of granular ma- 
terials. We do not see any evidence for hysteretic effects 
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